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ABSTRACT 

The two anomalous X-ray pulsars (AXPs) with well-sampled timing histories, IE 
1048.1—5937 and IE 2259+586, are known to spin down irregularly, with 'bumps' 
superimposed on an overall linear trend. Here we show that if AXPs are non-accreting 
magnetars, i.e. isolated neutron stars with surface mag netic fields Bq ^ 10^° T, then 
they spin down electromagnetically in exactly the manner observed, due to an effect 
called 'radiative precession'. Internal hydromagnetic stresses deform the star, creating 
a fractional difference e = (I3 — ^ 10^^ between the principal moments of inertia 

Ii and I3; the resulting Eulerian precession couples to an oscillating component of the 
electromagnetic torque associated with the near-zone radiation fields, and the star 
executes an anharmonic wobble with period Tpr ^ 27r/ef](t) ~ 10 yr, where Q(t) is the 
rotation frequency as a function of time t. We solve Euler's equations for a biaxial 
magnet rotating in vacuo; show that the computed ^}{t) matches the measured timing 
histories of IE 1048.1-5937 and IE 2259+586; predict for the next 20 years for 
both objects; predict a statistical relation between (dQ/dt) and Tpr, to be tested as the 
population of known AXPs grows; and hypothesize that radiative precession will be 
observed in future X-ray timing of soft gamma-ray repeaters (SGRs). 
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1. INTRODUCTION 

Anomalous X-ray pulsars (AXPs) are a sub-class of X-ray pulsars distinguished by pulse 
periods lying in a narrow range between 6 s and 12 s, low X-ray luminosities, soft X-ray spectra, 
no detected optical counterparts, no detected orbital Doppler shifts of pulse arrival times, and 
associations with shell-type supernova remnants (Mereghetti & Stella 1995; van Paradijs, Taam & 
van den Heuvel 1995). The nature of AXPs is uncertain; possibilities include an accreting neutron 
star in a binary, with a low-mass white dwarf or He-burning star as a companion (Mereghetti & 
Stella 1995; Baykal & Swank 1996; Mereghetti, Israel & Stella 1998), an isolated neutron star 
accreting from a residual disk following a phase of common-envelope evolution (van Paradijs et al. 
1995; Ghosh, Angelini & White 1997), a non-accreting massive white dwarf (Morini et al. 1988; 
Paczyiiski 1990; Usov 1994), and a non-accreting magnetar (Thompson & Duncan 1996; Heyl &: 
Hernquist 1998; Kouveliotou et al. 1998). 

The two AXPs with well-sampled timing histories, IE 1048.1-5937 and IE 2259-^586, are 
observed to spin down irregularly: the rotation frequency Q of each object decreases linearly 
with time t on average, but there are 'bumps' superimposed on this average trend during which 
= dQ /dt < fluctuates by a factor of between two and five every five to ten years (Mereghetti 
1995; Baykal & Swank 1996; Baykal et al. 1998; Oosterbroek et al. 1998; and references therein). 
In accreting-star models of AXPs, the average spin-down trend is attributed to the accretion 
torque acting on a near-equilibrium rotator, and the bumps are the result of white torque noise 
(e.g. due to disk inhomogeneities) as in ordinary binary X-ray pulsars (Baykal &: Swank 1996). 
In isolated-star models, the spin-down trend is attributed to magnetic-dipole braking, and the 
bumps are analogous to glitches observed in rotation-powered pulsars like Vela (Usov 1994; Heyl 
& Hernquist 1998). 

In this Letter, we present unequivocal new evidence that AXPs are non-accreting magnetars. 
In fQ, it is shown that as a magnetar spins down it wobbles anharmonically, with a period of five 
to ten years, due to an effect called 'radiative precession'. The spin-down signature of the wobble, 
calculated theoretically in §^ matches closely the bumpy timing histories of IE 1048.1—5937 and 
IE 2259-1-586. The theory is used to predict 0(i) over the next 20 years for both objects and 
yields a statistical relation between bump recurrence time and average spin-down rate for the AXP 
population as a whole. Implications for the internal and magnetospheric structures of magnetars, 
including AXPs and soft gamma-ray repeaters (SGRs), are explored in 

2. RADIATIVE PRECESSION 

A magnetar is a triaxial body in general. Hydromagnetic stresses arising from non-radial 
gradients of the superstrong internal magnetic field, e.g. between the magnetic poles and equator 
if the field is dipolar, deform the star, inducing matter-density perturbations Sp ~ Bf^/ ^lqc^ 
and hence a fractional difference e 6pR^/Ii ~ 2 x 10~^ (-6^/10^'^ T)^ between the principal 
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moments of inertia (Goldreich 1970; de Campli 1980; Melatos 1998). Here, Sin is the characteristic 
magnitude of the internal magnetic field, Cg is the isothermal sound speed (cg ~ 3~^/^c), and 
R is the stellar radius; one has B^^^ ~ Bq if the internal field is confined to the stellar crust 
and Sin ^ -Bq if it is generated deep inside the star, e.g. in a convective dynamo (Thompson Sz 
Duncan 1993). In a rotation-powered pulsar with Bq ^ 10^ T, the hydromagnetic deformation is 
comparable to the elastic deformation arising from shear stresses in the crystalline stellar crust, 
and the principal axes of inertia are oriented arbitrarily with respect to the axis m of the external 
magnetic dipolc (Goldreich 1970; dc Campli 1980; Melatos 1998). In a magnetar, however, the 
hydromagnetic deformation is much larger, and m is approximately parallel to one of the principal 
axes (e3, say); the alignment is not exact due to the complicated structure of the internal field 
near its generation site, cf. the non-axisymmetric, magnetically modified Taylor columns seen in 
simulations of the geodynamo (Glatzmaier & Roberts 1996). Provided that the rotation axis ft 
is not parallel to m ~ 63, as is usually the case for rotation-powered pulsars, the star precesses 
about 63 with period Tpr - 27r/eJ^ 85(5in/10^° T)-2(J^/lrads-^)-i yr. 

The Eulerian precession is not free. It couples to a component of the vacuum magnetic-dipole 

torque associated with the asymmetric inertia of the near-zone radiation fields outside the 
magnetar; the electromagnetic energy density (and hence inertia) of the near-zone radiation 
fields is greater at the magnetic poles than at the equator by an amount ^ Bq/hq, resulting 
in an oscillatory, precession-inducing torque (Goldreich 1970; Good & Ng 1985; Melatos 1998). 
The near-field torque exceeds the familiar braking torque (oc J7^) by a factor c/flR S> 1 and 
acts on a commcnsuratcly shorter time-scale, Tnf tq^R/c ~ 6(i?o/10"^'^ T)^^(r2/1 rads~^)~^ yr, 
where tq = ^iqc^Ii/2'kBIR^VL^ « 2 x 10^(So/10^° T)-2(17/i rads~i)"Vr is the characteristic 
electromagnetic braking time. Since the near-field torque is directed along O x m, it does not 
change Q. for a spherical star. Nor does it change O for an aspherical star, provided e is large 
enough to give Tpr ^ Tnf, so that the near-field torque averages to zero over one precession period. 
When the dominant deformation is hydromagnetic, however, one finds Tpr /Tnf ~ 14(Bo/-Bin)^, close 
to unity (i.e. strong coupling) provided i?in is moderately larger than Bq as expected (Thompson 
& Duncan 1993). Under these circumstances, the star executes an anharmonic wobble, called 
radiative precession, with period Tpr (« Tnf), and the near-field torque changes on the precession 
time-scale (Melatos 1998). 



3. THEORY OF BUMPY SPIN-DOWN 

3.1. Solution of Euler's Equations of Motion: Past and Future 0(t) 

We now show that the timing signature of radiative precession matches the observed bumpy 
spin-down of AXPs by solving numerically Euler's equations of motion for a biaxial, dipole magnet 
rotating in vacuo, 

Cli = — eri2^^3 + ^"^(T^ cosx -I- sinx) -I- 602(^^1 sinx -I- f23COSx)], (1) 
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^2 = eilir^s + ^0 ^''"o ^["'^^^^2 + ^1 cosx + sinx)(^^i sinx + ^^3 cosx)]; (2) 
= — JIq ^Tq"^ sinx[af^^(— J^i cosx + f^s sinx) + &f^2(^i sinx + cosx)]- (3) 

In (|l|)-(^), subscripts denote components along the principal axes of inertia, x is the angle between 
m and e^, and we have a = 0.33, b = 0.094c/r2oi?, and 0,q = 0,{t = to), where to is an arbitrary 
origin; for derivations of the equations, see Goldreich (1970) and Melatos (1998). Terms in (@)~(^) 
proportional to e give rise to Eulerian precession, terms proportional to bTQ^ are associated with 
the near-field torque, and terms proportional to otq"^ cause secular braking. A biaxial, dipole 
magnet rotating in vacuo is an idealized model of a hydromagnetically deformed magnetar. In 
reality, such a body is triaxial (if it is indeed rigid), has high-order and/or off-centered multipoles 
contributing to the near-zone magnetic field, and is surrounded by a plasma magnetosphere. The 
values of a and b refiect, in a coarse way, the magnetization state of the stellar interior and the 
distribution of magnetospheric currents (Melatos 1998). 

In Figures Q and |2|, we plot the computed r2(t) from (|ll)~(|3|) on top of the X-ray timing 
histories of IE 1048.1—5937 and IE 2259+586 respectively, extending the theoretical curves 20 
years beyond the present as a testable prediction. The timing history of IE 2259+586 is sufficiently 
well-sampled that only one good fit is possible. In the case of IE 1048.1—5937, the solid curve is 
the favored fit, with Tpr ~ 18 yr, but an alternative fit, with similar average slope and Tpr ~ 6yr, is 
also acceptable. Multiple good fits are hard to find. 

Euler's equations of motion contain just three unknown parameters: e, tq, and x- It is 
significant that the theory agrees with the observations as well as it does despite its idealized 
nature — particularly as only two of the three parameters are truly free, since one needs 
strong torque coupling (rpr ~ t^u or equivalently eiloTQ ~ c/QqR) in order to get any bumps 
at all. Moreover, the best agreement with observations is achieved for parameter values that 
are consistent with general physical considerations. If AXPs are hydromagnetically deformed 
magnetars with Bi^ Bq ^ 3 x 10^^ T, one expects e 10~^, OqTo ~ 10^^, and x relatively small 
(cf. X ~ 11° for the Earth), as discussed above. 

The theoretical curves do not pass exactly through every available data point, and more 
departures are expected in the future, e.g. the slope between the last two data points in Figure || 
is ~ 0.6 times the solid-curve theoretical slope at that epoch. Indeed, a formal estimate of the 
chi-square for the solid curve in Figure ||, taking published timing uncertainties at face value, 
implies that the fit is poor: one finds a chi-square of ~ 8 x 10^ with 10 degrees of freedom, notably 
inferior to multiple-glitch models, for example (Heyl & Hernquist 1998). This is because a biaxial, 
dipole magnet is an over-idealized model of an AXP, as noted above; the chi-square likelihood 
improves dramatically for a more realistic model with just two extra parameters, e.g. a non-dipolar 
near-zone magnetic field and a triaxial ellipsoid of inertia. However, our aim in this paper is not 
to model Q(t) in detail, but to account for key gross features of the data — the average spin-down 
rate, bump recurrence time, and bump amplitude — with a simple physical theory. In this regard, 
the agreement with observations is good. 
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One data point in the timing history of IE 2259+586, at t — to = 15.4 yr, represents a spin-up 
event lasting at most 0.8 yr at the la level of uncertainty (Baykal & Swank 1996). Spin-up cannot 
be explained by radiative precession because (i) (l is always negative, even while decreases by 
up to a factor of five during a bump, and (ii) there is no natural way to explain a 0.8 yr time-scale. 
If taken at face value, the spin-up event must be a different phenomenon, e.g. a discontinuous 
change of internal structure analogous to the glitches of rotation-powered pulsars like Vela (Usov 
1994; Heyl & Hernquist 1998). 

3.2. Predicted Population Statistics 

Numerical studies show that the bump recurrence time and average spin-down rate satisfy 
Tpr oc B^^Q~^ f{x) and (17) oc BQ^l'^g{x), with 1 < /(x), g{x) < 10. In other words, the greater 
the magnetic field of an AXP, the shorter is its precession period and the greater is its spin-down 
rate, with 

{i}) ^-2x 10~^(5o/-Bin)^(f7/lrads"^)2(rpr/lyr)~irads~Vr"^- (4) 
In addition, one finds that Q increases from ~ (0,) to ~ over a time ~ 0.25rpr during the 
course of a bump, yielding a bump amplitude 

AQpr/^ ~ 5 X 10--'{Bo/Bir,f{n/lrads-^) ; (5) 

this number is similar for all AXPs. Both the relations (^) and (|5|) will be subject to observational 
testing as the population of AXPs with measured timing histories swells over time. Note that 
they are statistical relations, with scatter expected about an overall trend. The detailed structure 
of the magnetic field inside an AXP — which does not enter into the idealized theory presented 
here, except through x — is likely to differ from object to object, affecting Tpr, (Cl), and AQpj. 
significantly, as the broad ranges of f{x) and g{x) attest. 

4. DISCUSSION 

Why is radiative precession not observed in rotation-powered pulsars with Bq ^ 10^ T? 
These objects are deformed hydromagnetically like magnetars, with an added elastic deformation, 
and they spin down electromagnetically. Young pulsars like the Crab, with tq ~ 10^ yr and 
nn/c ^ 10-^, ought to precess with period Tpr ~ 10 yr, yet there is no clean evidence of bumpy 
spin-down in radio timing data, nor of concomitant changes in pulse profile (e.g. relative height 
and separation of conal components) and polarization characteristics (e.g. position-angle swing). 
Lyne, Pritchard & Smith (1988) reported a quasi-sinusoidal variation in Crab timing residuals 
with a period of « 20 months but judged it likely to be an artifact of an overlooked glitch. 
(See also Melatos 1998.) The only reliable detection of pulsar precession to date has been the 
general-relativistic geodetic precession of the binary pulsar PSR B1913+16 (Weisberg, Romani &; 
Taylor 1989; Kramer 1998). 
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One possible explanation is that radiative precession is damped in pulsars with Bq ^ 10^ T. 
Frictional dissipation inside the star, due to time-dependent elastic strains in the crust and/or 
imperfect crust-core coupling, is thought to occur on a time-scale ^ lyr (Link, Epstein & Baym 
1993), rapidly aligning ii with 63. (In the case of the Earth, dissipation restricts fl to fluctuating 
within just 1" of 63 under the action of solar and lunar tides; see Bursa &; Pec 1993). In a 
magnetar, where the magnetic energy exceeds the mechanical energy of rotation, the stiffening 
effect of the superstrong magnetic field may hinder the elastic strains and/or sheared fluid flows 
responsible for internal damping. A second possible explanation is that conduction currents 
in the magnetosphere of a pulsar with Bq ^ 10^ T, where electron-positron pair production is 
plentiful, modify the electromagnetic torque in such a way as to suppress the precession-inducing 
near-field component. In a magnetar, where it is thought that pair production is quenched, e.g. 
by positronium formation (Usov Sz Melrose 1996; J. Heyl 1999, private communication) or photon 
splitting (Baring & Harding 1998), the vacuum magnetic-dipole torque, with its unmodified 
near-field component, is a closer approximation to reality. Both explanations imply that radiative 
precession will not be suppressed in SGRs, where one has Bo ~ 10^^ T (Kouveliotou et al. 1998). 

What does radiative precession teach us about strongly magnetized neutron stars themselves? 
Firstly, the close agreement between theory and observation implies that AXPs are indeed 
non-accreting magnetars, and that accretion is not needed to explain fluctuations in Q, contrary 
to claims in the literature (e.g. Mereghetti 1995). Secondly, if SGRs are magnetars too, as 
indicated by recent X-ray timing data (Kouveliotou et al. 1998), they ought to exhibit bumpy 
spin-down just like AXPs, perhaps punctuated by brief, glitch-like spin-up events during the 
gamma-ray bursts themselves. This constitutes a direct test of the magnetar model for SGRs. 
Thirdly, the fact that bumpy spin-down is seen in AXPs implies that B[^i is at most a few times 
Bq to ensure strong coupling, i.e. Tpr ~ Tnf, as discussed above. This is indirect evidence that 
the magnetic field of these objects is generated relatively near the stellar surface. Fourthly, 
radiative precession places an upper limit on the fluid viscosity r/ of a newly born magnetar; if rj 
is too high, ri aligns with 63 ~ m before the stellar crust crystallizes, and there is no precession 
(Melatos 1998). The upper limit obtained in this way — that the viscous damping time is greater 
than the crust crystallization time of ~ 1 yr validates certain semi-quantitative calculations 
of r] for pulsars with Bq ^ 10^ T by Cutler & Lindblom (1987). We remark in closing that the 
predicted hydromagnetic deformation (e ~ 10~^) constitutes a misaligned mass quadrupole which 
generates gravitational radiation. However, the radiation is too weak to be detected by planned 
gravitational- wave interferometers like LIGO and VIRGO, because AXPs are slow rotators and 
the gravitational- wave amplitude is a strongly rising function of fl. 

I am indebted to Lars Bildsten for bringing the existence of spin-down irregularities in AXPs 
to my attention, and for subsequent discussions. I also thank Jon Arons, Peter Goldreich, and the 
referee, Jeremy Heyl, for comments. This work was supported by NASA Grant NAG5-3073, NSF 
Grant AST-95-28271, and by the Miller Institute for Basic Research in Science through a Miller 
Fellowship. 
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Fig. 1.— Rotation frequency Q versus time t for the AXP IE 1048.1-5937, with to = JD 2444000. 
The squares and accompanying la error bars are X-ray timing measurements made over the last 
20 years by the satelhtes Einstein, EXOSAT, Ginga, ROSAT, ASCA, RXTE, and BeppoSAX 
(Oosterbroek et al. 1998, and references therein). The sohd curve is the solution to Euler's equations 
of motion (|)-(|) for = 0.976 rads-^ I^qTo = 3.8 x 10^°, e = 6.4 x 10-^ and x = 3.5°, with 
initial conditions Jli^o = 0.476ilo and 1^2,0 = — 0.568$7o- The broken curve is the solution for 
Qq = 0.9762 rads~^ QoTq = 3.5 x 10^°, e = 6.3 x 10^^, and x = 18°, with the same initial 
conditions. The initial conditions determine the initial phase of the oscillation and are otherwise 
insignificant. 
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Fig. 2.— Rotation frequency n versus time t for the AXP IE 2259+586, with to = JD 2443000. 
The squares and la error bars are X-ray timing data from the satehites cited in Figure |l] as weU 
as from HE AO 1 and Tenma (Baykal et al. 1998, and references therein). The sohd curve is the 
solution to Euler's equations of motion (|l|)-(|) for JIq = 0.900356 rads~\ QqTo = 2.35 x 10^^, 
e = 3.4 X 10~®, and x = 13°, with initial conditions 1^2,0 = — 0.668r2o and Q^fi = 0.658r2o- 



